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ABSTRACT 


Effects of surface heat transfer on boundary-layer transition 
are analyzed in a three-part study using the Saffman-Wilcox 
transition model. In the first part of the study, model pre¬ 
dictions are compared with experimental data for cooled and 
heated aerodynamic boundary layers on smooth flat surfaces and 
for cooled aerodynamic boundary layers near the stagnation 
point of a roughened blunt body. Consistent with measure¬ 
ments, tne model predicts, on the one hand, that heating 
destabilizes a smooth-surface aerodynamic boundary layer 
and, on the other hand, that cooling destabilizes a rough- 
surface aerodynamic boundary layer. Differences between 
predicted and measured transition-point locations are with¬ 
in experimental error bounds. Then, incipient transition 
conditions arc determined for a small, heated hydrodynamic 
body. Again model predictions agree with measurements which 
indicate that relatively small amounts of surface heating 
have a strong stabilizing effect on hydrodynamic boundary 
layers. In the final part of the study, transition location 
is determined for a large hydrodynamic body; results indi¬ 
cate that large surface heating rates are not substantially 
more effective than smaller rates., i 
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NOTATION 


SYMBOL 


DEFINITION 


e 

f r 

F(n) 
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00 
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U e 
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Re 


Pressure coefficient, (p-p^J/C^/^pU^ 2 ) 

Specific heat at constant pressure 
Specific turbulent energy 

°arsons-Goodson—Goldschmied shape parameter 
Nondimensional self-similar velocity profile 
0 for two-dimensional flow; 1 for axisymmetric flow 
Roughness height 

Parsons-Goodson-Goldschmied shape parameter 
Pressure 

Laminar and turbulent Prandtl numbers 
Surface heating rate; average value of q w 
Radial distance from symmetry axis 
Parson-Goodson-Goldschmied shape parameters 
Empirical constant 

Reynolds number based on arclength, plate length 

Reynolds number based on momentum thickness 
Incipient transition Reynolds number based on V 
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Parsons-Goodson-Goldschmied shape parameter 
Parsons-Goodson-Goldschmied shape parameter 
Temperature 

Freestream total temperature 
Freestream turbulence intensity 
Velocity component in x,y direction 
Boundary-layer-edge velocity 
Freestream velocity 
Body volume 
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P 
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Distance parallel to, normal to body surface 

Parsons-Goodson-GoIdschmied shape parameter 

Empirical parameters 

Empirical constants 

Similarity variable 

Momentum thickness 

Karman's constant; thermal conductivity 
Empirical parameters 

Modified Polhausen pressure gradient parameter 

Molecular viscosity 

Kinematic viscosity 

Fluid density 

Empirical constants 

Specific dissipation rate 


subscripts 
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Surface 
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1. INTRODUCTION 


In recent years, substantial progress has been made toward 
maintaining laminar flow, and hence low drag, over aerodynamic/ 
hydrodynamic bodies at practical flow speeds. In the case 
of aerodynamic bodies, the use of carefully designed air¬ 
foil shapes (e.g. to provide strong favorable pressure gra¬ 
dient) and the use of boundary-layer-control techniques 
(e.g. suction) has successfully delayed transition to chord- 
length Reynolds numbers, Re St , of the order of 10 million. 

Even larger values of Re s ^_ have been achieved for small hydro- 
dynamic bodies through the use of surface heating and pres¬ 
sure gradient. For hydrodynamic bodies, extrapolations based 
on linear stability theory indicate that, with practicable 
amounts of surface heating, values of Re St in excess of 200 
million may be possible on relatively large hydrodynamic 
bodies 1 . 

The reduction in drag which can be achieved by maintaining 
laminar flow over any vehicle is attractive because of the 
reduced power requirements to move the vehicle. However, a 
penalty is generally paid in maintaining laminar flow on 
lifting bodies in that a laminar boundary layer separates 
much more easily than does a turbulent boundary layer, re¬ 
sulting in a significant reduction in lift. Hence, main¬ 
taining laminar flow is practical mainly for nonlifting 
bodies. Submarines and torpedos or, more generally, hydro- 
dynamic bodies fall into the latter class. 

This study focuses on the observed pronounced effects of 
surface heat transfer on boundary layer transition. Most 
importantly, this report includes transition predictions 
based on a relatively new transition theory for small and 
large heated hydrodynamic bodies. 
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Section 2 presents the transition equations, including a 
modification needed to improve transition-prediction ac¬ 
curacy when surface heat transfer is present. Included in 
Section 3 are transition computations for various aero¬ 
dynamic and hydrodynamic boundary layers. The concluding 
section summarizes results and conclusions. 





















2. FORMULATION 


2.1 EQUATIONS OF MOTION 

The Saffman-Wilcox transition model 2 ^ is the basic tool 
used in this study to analyze effects of surface heat 
transfer on boundary-layer transition. The model's accuracy 
has previously been demonstrated for a wide variety of flows 
ranging from incompressible boundary layers to hypersonic 
blunt-body flows. These applications have tested the model's 
ability to predict transition sensitivity to effects of free- 
stream turbulence, suction, surface roughness and pressure 
gradient. In all cases, accuracy acceptable for most en¬ 
gineering purposes has been obtained. 

For incompressible boundary layers (l.e., for very small Mach 
numbers), the equations of motion which constitute the tran- 

sition model are: 

Mass Conservation 


H <rJ u> + W (r,Iv> 


(1) 


Momentum Conservation 


3u . 3u _ 

o u 37 + pv 37 ‘ 


d£ , 1 JL [*rJ (u + -) —1 

' -dS + J 9y L n 9y J 


( 2 ) 


Energy Conservation 
3(C T) 3(C T) 

+ pv -5r" 
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Turbulent Energy 


3e + 

pu ^ + 


pv 


0 = [e.*l0|-6*Pfl]pe + ^^[ r3(lJ + a *l ) 3y] 


(U) 


Turbulent Dissipation Rate 


pu 


II 1 + pv = [cxlf^l - 6pn]p£2 2 - + ] 


(5) 


m Equations 1-5, x and y are orthogonal coordinates parallel 
to and normal to a body surface. The quantity r is e ra 
coordinate from the body's symmetry axis while .1=0 for wo 
dimensional flow and J=1 for axisymmetric flow. The veloc 
components In the x and y directions are denoted by u an v, 
fluid density, temperature, pressure, specific heat an, 
viscosity are denoted by p, T, p, C p , and u respec 
quantities e and 0 are specific turbulent energy and .peomc 
turbulent dissipation rate; their ratio, e/Sl, is e ur 
eddy viscosity. 

Both air and water are considered in this study. The various 
thermodynamic properties for air are related through the per¬ 
fect gas law and the Sutherland viscosity law. ppen 
lists the pertinent thermodynamic properties of water. 

Several empirical parameters appear in Equations 1-5. For 
transitional flows, past studies 2 ' 4 have established the 

following values: 


Pr T * 0.89 


6 = 0.15 
a = 0 .50 » 


, B* = 0.09 
a* = 0.50 


(6) 






(7) 


a* = 0.30 [1 - (1 - X*) (1 - Re T /R 0 ) *' " Re T /R o ) -' 

a= .1638 [1 - (1 -X)(l - Re T /R Q )H(l - R e T /R Q )] (8) 

where H(x) is the Heaviside stepfunction and Re^ is turbulent 
Reynolds number defined by 

Re^ = e/fiy (9) 

i 

Finally, the parameters A# and R q are constants whose values 
are 


a * = 0.105 , : > 0 = 0.10 


( 10 ) 


The quantity A depends upon the freestream turbulence level, 
T*, and pressure gradient parameter, A, defined by 


T' 



and 


^e dj_ dU. 
P w v e dx 


( 11 ) 


( 12 ) 


The functional dependence of A upon A and T' is 

A= . I05(l + 2H(A) j 1 - exp [--40 A ] } exp [-3T' 2 ]) ^3) 

Equations 1-3 are the time-averaged conservation equations 
with classical eddy-viscosity and eddy-heat-diffusivity clo¬ 
sure approximations. Equations 4 and 5 are -onlinear diffusion 
equations for e and JJ which provide a description of the growth 
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of disturbances in a laminar flow up to, through, and beyond 

2 

transition. Wilcox presents a thorough discussion of the way 

in which the model equations are used to predict boundary- 

layer transition. Boundary conditions suitable for flow over 

both smooth and rough surfaces are given by Wilcox and 
h 

Chambers . 


2.2 NEUTRAL STABILITY CONSIDERATIONS 


Qualitative features of model-predicted transition can be 
conveniently determined by dropping the convection and dif¬ 
fusion terms in the turbulent-energy equation (Equation 4). 
In doing this we can determine the neutral-stability point 
which is defined as the point in a boundary layer where tur¬ 
bulent energy generation, ot#pe | 3u/9y | , just balances turbu¬ 
lent energy dissipation, 8*p 2 fte. Hence, neutral stability 
is defined by the following condition: 


m?x 


a* 

B* 


9u /3y j n 1 

pft 1 


( 1*0 


Then, noting for laminar boundary layers that ft and a* are 
approximately 


ft 


20 y 
6 P^y 2 


(15) 


ct#= 0.30 A# 


Equation 14 simplifies to 


m f X 


y 2 |9u| _ 40 
v 197 ! ~ ITT 


(lo) 


(17) 
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For zero-pressure-gradient boundary layers in air, the velocity 
distribution can be written as 


u = U e F(n) (18) 



Using Equations 18 and 19, Equation 17 simplifies to the 
following: 

( 20 ) 

Now for incompressible flow, F(n) is only weakly dependent 
upon surface heat transfer. Hence, Equation 20 implies that 
the neutral stability Reyi olds number, U e x/v w , varies inversely 
as A* 2 . Experiments 5 indicate that transition Reynolds number, 

p 

U x,/v , is inversely proportional to (T /T ) . Assuming that 
e u w w tr 

the transition Reynolds number is proportional to the neutral 
stability Reynolds number (as it often is) implies that A* 
should be proportional to T w /T p . 

Computations with A* independent of T w /T £ verified that such a 
dependence is needed. The unmodified model predicted that 
heating (cooling) has a stabilizing (destabilizing) effect on 
aerodynamic boundary layers, in contrast to the measured de¬ 
stabilizing (.stabilizing) effect. These predictions were un- 
surprising as Shamroth and McDonald^ find a similar reversal 
in predicted transition sensitivity to heating and cooling 
with their turbulence-model transition method. Shamroth and 
McDonald resolve the problem by making a parameter similar to 









X* an increasing function of T/T . While this model re- 

W 0 

vision is sufficient to yield accurate transition location 
for aerodynamic boundary layers, it is inadequate for the 

hydrodynamic case. That is, according to Equation 20, tran- 

2 

sition location, x. , is proportional to y /X# . For a 

b _ 0 W 

liquid, y varies as T “ near room temperature (see Appen- 

dix «. ..ence, assuming X* is proportional to T.. would ir.rlv 

that x varies as T Therefore, increasing surface ten- 

t w 

perature would decrease , i.e., heating would destabilize 
a hydrodynamic boundary layer. Since heating stabilizes 
such a boundary layer, a different modification is clearly 
needed. 

The parameter X* should more appropriately depend upon y^/y^. 
On the one hand, y /y increases as T /T increases for air, 

W G W 0 

while, on the other hand, y, r /y Q decreases as T /T increases 
for water. Thus, if X* were proportional to some power of 
y w //y e’ the mode - 1 - cou l d accurately predict effects of sur¬ 
face heat transfer on both aerodynamic and hydrodynamic 

o/h 

boundary layer transition. Since y^T for air, the fol¬ 
lowing revised form of Equation .0 is proposed: 

x* = 0.105 (y w /y p ) i4/3 (10a) 

An additional, less obvious, modification is needed. Note 
that the parameter X partially control® the rate at which 
disturbances are amplified beyond the neutral stability 
point. Furthermore, if X ever becomes sufficiently large 
(relative to X*) so that dissipation-rate production, 
apft 2 |9u/9y|, overtakes dissipation-rate dissipation, Bp 2 ^ 3 , 
(see Equation 5) before the neutral stability point is 
reached, transition may never occur. Thus, if X remains 
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unaltered, the possibility exists that a finite amount of 
cooling (heating) will cause an aerodynamic (hydrodynamic) 
boundary layer to forever remain laminar, a physically un¬ 
realistic prediction. Hence, Equation 13 must be repxaced 
by: 

X=. 105(y w /y p ) i</3 { 1 + 2H(A)[1 - exp(-40A) lexp(-3T’ 2 )} (13a) 

All of the computations presented in the next section hav^ 
been performed using the transition model defined by Equations 
1 through 9, 10a, 11, 12, and 13a; as noted earlier, appro¬ 
priate surface boundary conditions are given by Wilcox and 

Chambers 1 *. 
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3. APPLICATIONS 


In the first part of this section, the transition model is 
used to predict transition sensitivity to surface heat trans¬ 
fer for smooth- and rough-surface aerodynamic boundary layers; 
computed transition-point locations are compared with cor¬ 
responding experimental data. Then, effects of surface 
heating on a small hydrodynamic body are computed; qualita¬ 
tive comparisons are made with experimental data. Finally, 
transition location is predicted on a large hydrodynamic 
body. 

HEATED AND COOLED AERODYNAMIC BOUNDARY LAYER^. 

One of the easiest of all flows to analyze is the incompres¬ 
sible flat-plate boundary layer (FPBL). Furthermore, analyz¬ 
ing this flow provides a good test of the transition model 
as detailed measurements have been made to determine transi¬ 
tion sensitivity to surface heat transfer. Using an incompres¬ 
sible version of DCW Industries' EDDYBL computer code , tran¬ 
sition computations were performed for an incompressible PPBL 

with 

0.5 < T/T < 3.0 (21) 

As in all computations in this study, the value of ft at the 
boundary-layer edge, , was given by 

= .Ol85U*/y e ( 22) 

a value generally used in EDDYBL transition calculations. Re¬ 
sults of the computations are shown in Figure 1; experimental 
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data of Zysina - Molozhen and Kuznetsova are included in 
the figure for comparison. 

Because the freestream turbulence level in the experiments was 
not available, computations were done for a relatively high 
level (T'=l.25$) and for a relatively low level (T'=.03$). As 

shown in the figure, freestream turbulence level has only a 
slight effect on transition Reynolds number based on surface 
conditions, (Re Xt ^ w sU e x t //v w‘ Consistent with the measure¬ 
ments, computed (Re Xt ) w is approximately inversely propor¬ 
tional to (T w /T e ) 2 ; computed values of (Re xt ) w generally are 
within about 20$ of corresponding measured values. 

A second, more subtle, test of the theory is to apply the 
equations to rough-wall aerodynamic ooundary layers. Surface 

cooling can actually reverse its stabilizing role to one of 
destabilization when the surface is rough. The explanation 
for this phenomenon is well understood. The destabil zation 
occurs because cooling thins the boundary layer, thus making 
the roughness look larger relative to a boundary-layer thick¬ 
ness such as momentum thickness, 6. Hence, the surface looks 
rougher, whereby transition occurs earlier. 

To test the model's ability to predict this phenomenon, com¬ 
putations were performed for Mach 5 airflow past a roughened 
spherical body with nose radius of 2.5 inches; roughness 
heights, k, between 1.5 mils and 10 mils were used in the 
computations. Freestream unit Reynolds number ranged from 
3 to 10 million per foot. (The compressible version of the 
model equations was used for these calculations). The ratio 
of surface temperature to freestream total temperature, T. , 
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was varied from 0.2 to 0.8. Figure 2 compares results of 

the computations with a correlation of experimental data'. 

Consistent with the measured effect, momentum thickness 

Reynolds number at transition, Re Q , increases with de- 

°t 

creasing surface temperature. 



5 10 20 



Figure 2. Comparison between computed and measured 
effects of surface cooling on rough-wall 
aerodynamic boundary-layer transition. 
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HEATED HYDRODYNAMIC BOUNDARY LAYERS ON SMALL BODIES 


The second round of applications is to a small heated hydro- 
dynamic body. The body, shown in Figure 3, is a Parsons- 
Goodson-Goldschmied 8 (PPG) minimum-drag hull shape which will 
be referred to as the H-2 body; the PPG body parameters are 
given in Appendix B. The objective of the computations was 
to determine incipient transition conditions for various 
surface temperatures. By definition, incipient; transition 
occurs when the transition point is located at the maximum 
body radius. 


All computations were performed with the following conditions 
specified: 


Ambient Temperature, T 0 
Roughness Height, k 
Freestream Turbulence, 
Freestream Dissipation Rate, G e 


= 71°F 
= 16 yin 
= . 01 % 

= .0185 Ug/y g 


The pressure distribution was obtained from the Douglas- 
Neumann 9 potential flow program; Figure 4 shows the computed 
pressure coefficient, Cp, defined by 


■p = (p-p„)/(i/2pu£> 


(23) 


In the figure, z is axial distance from the stagnation point. 
Three series of computations were performed in which surface 
temperature was held constant at values of 71 f> 7*3 F> and 
8 l°F. Incipient transition Reynolds number, Re^. (based on 
freestream flow conditions and V 1/3 where V is body volume', 
was computed for each value of T w> Results of the calculations 
are summarized in Table 1; in addition to Re y , the table gives 
freestream velocity, U w , transition Reynolds number based on 
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Figure Computed pressure coefficient for 

the small hydrodynamic body 





In Equation 24, s = 4 feet (the maximum radius point), and 

1 / 

q (s) is the local heat transfer. Figure 5 shows computed Re v 
w * 

as a function of T -T and q . As can be seen from the figure, 

we w 

heating stabilizes the boundary layer and hence increases the 
value of Rey This prediction is qualitatively consistent 
with the observed stabilizing effect of surface heating when 
the fluid is a liquid. The magnitude of the effect is also 
realistic" 1 ". 


Table 1. Computed Incipient Transition Conditions for 
the H-2 Body for Varying Surface Temperature. 


- ,kwatts^ 
q w ( ' -ft 2-) 



U (knots) 

oo 

10 _6 Re v 

25.3 

9.4 

35.0 

13.0 

47.4 

17.6 


10 -6 Re, 


18.2 
24.5 
33.8 


An interesting feature of the predictions is revealed in 

Figure 5- Specifically, Re v increases approximately linearly 

with T w ~T e for tne range of temperatures considered. However, 

Re,, increases much less rapidly with q . Hence, increasing 
V w 2 

the heating rate beyond about 2 kwatts/ft*" may not yield 

substantial increases in Rey 
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3.3 HEATED HYDRODYNAMIC BOUNDARY LAYERS ON LARGE BODIES 

The final application of the transition model is to a large 
heated hydrodynamic body. Again, a PPG hull shape, referred 
to as .he R-9a body, is used (see Appendix B). The objective 
of these computations was to determine transition point lo¬ 
cation for specified flow conditions. 

Two computations were performed; in one, the freestream flow 
speed, U B , was 35 knots and in the other U OT was 30 knots. In 
both computations, the following conditions were imposed: 

Ambient Temperature, T g 
Surface Temperature, T w 
Roughness Height, k 
Freestream Turbulence, T' 

Freestream Dissipation Rate, 

Figure 6 shows the Douglas-Neumann-computed pressure distribu 

tion. 


* 55°F 
= 85°F 
= 16 yin 
= . 01 % 

= .0185 U 2 /u e 


Table 2 and Figure 7 summarize results of the computations. 
For both flow speeds, transiton occurs well upstream of the 
maximum radius. Computed arclength transition Reynolds num¬ 
bers, Re s , are not much larger than the value achieved on 
the H-2 body with q w = 1.30 kwatts/ft 2 (see Table 1). Noting 
the heating rates involved, the predictions are consistent 
with results of the preceding subsection. 
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4. SUMMARY AND CONCLUSIONS 


Results reported in Section 3 demonstrate the transition 
model’s accuracy for boundary layers with surface heat 
transfer, provided nonuniform viscosity effects are taken 
into account (Subsection 2.2). With no parameter adjustment, 
the model accurately predicts transition Reynolds number 
for (a) incompressible smooth-wall aerodynamic boundary 
layers with surface heating and cooling, (b) cooled rough- 
wall boundary layers on blunt bodies in a hypersonic air- 
stream, and (c) heated hydrodynamic boundary layers on small 

bodies. 

Results of the hydrodynamic computations agree with the 
measured strong boundary-layer stabilization attending 
small amounts of surface heating. The model also predicts 
that large amounts of heating are not significantly more 
effective in delaying transition than the smaller rates 
considered in the H-2 body computations of Subsection 3.2. 
The larger R-9a body computations support this prediction. 


. 







APPENDIX A 

THERMODYNAMIC PROPERTIES OF WATER 

The following thermodynamic properties of water are perti¬ 
nent to boundary layer transition: mass density (p;, speci¬ 
fic heat (C ), thermal conductivity, (<), molecular viscosity 
(y), rnd laminar Prandtl number (Pr L ). Values used in this 
study for these quantities have been obtained from . chlicht- 
ing 10 and are valid for temperatures ranging from about 40°F. 
to 110°F. 

Mass density, specific heat, and thermal conductivity are 
approximately constant and have the following values: 


= 1.936 

lbf•sec 2 /ft 4 

(Al) 

= 1 

Btu/lbf•°F 

(A2) 

= 0.35 

Btu/ft•hr•°F 

(A3) 


Viscosity and Prandtl number, by contrast, are strongly tem¬ 
perature dependent over this range of temperatures. An ap¬ 
proximate polynomial fit to the Schlichting data was used to 
calculate y in the computations of Section 3, i.e., 

-6 - 8 

y = 7.943*10 (T/600) lbfsec/ft 2 (A4) 

with temperature, T, in degrees Rankine. Finally, the laminar 
Prandtl number follows from its definition, namely, Pr L =yC p /<, 
wherefore, 

Pr T = 2.5KT/600)" 8 (A5) 

U 

Figure A1 compares Equations A4 and A5 with the measured de¬ 
pendences of y and Pr^ upon temperature. 


* 
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APPENDIX B 

HYDRODYNAMIC BODY SHAPE PARAMETERS 


The two hydrodynamic bodies analyzed in this study are 
specified in terms of the eight dimensionless parameters 
used by Parsons, Goodson, and ^cldschmied . The values of 
these parameters are as follows: 

Parameter H-2 Body R-9a Body 


f 

4.450 

5.000 

r 

x m 

0.427 

0.600 

k. 

1 

2.715 

4.000 

r n 

0.260 

0.050 

r i 

0.740 

0.700 

s i 

1.650 

1.750 

X i 

0.595 

0.800 

t 

0.180 

0.150 


Additionally, the H-2 body is 9-457 feet long and has a volume 
of 14.5 cubic feet; the R-9a body is 55 feet long and has a 
volume of 2500 cubic feet. 
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